CLAIMS 




1. (amended) A method for increasing the speed of line-by-line calculations of multiple 
overlapping absorption lines, said method comprising: 

(1) a step of dividing a domain of a Voigt function representing the shape of an 
absorption line into a first range around the peak of the Voigt function and a skirt 
portion not contained in the first range, replacing the first range with a cubic function 
whose function values and derivative values match with those of said Voigt function at 
the connection points with said Voigt function, and calculating the values and 
derivative values of said cubic function and the Voigt function in the skirt portion for 
each of first predetermined intervals; 

(2) a step of adding together the results of step (1) for a plurality of absorption lines; 

(3) a step of interpolating the results of step (2) to calculate the function values and 
derivative values in second predetermined intervals smaller than said first 
predetermined intervals; 

(4) a step of dividing said first range into a second range near the peak and a skirt 
portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
second cubic function whose function values and derivative values match with those 
of said "function representing the difference between the Voigt function and said cubic 
function" at the connection points with said "function representing the difference 
between the Voigt function and said cubic function", and calculating the values and 
derivative values of said second cubic function replacing said second range and the 
"function representing the difference between the Voigt function and said second 
cubic function" in the skirt portion for said second predetermined intervals; 
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(5) a step of adding the results of step (4) to the results of step (3) for a plurality of 
absorption lines; 

(6) a step of interpolating the results of step (5) to calculate the function values and 
derivative values for intervals smaller than said second predetermined intervals; and 

(7) a step of adding the values of the "function representing the difference between 
the Voigt function and said second cubic function" to the results of step (6) for a 
plurality of absorption lines in said second range. 

2. (amended) A method in accordcince with claim 1, wherein (8) the step of calculating 
function values and derivative values by interpolation is repeated while narrowing the 
interval until it becomes an interval of a minimum unit. 

3. (amended) A method in accordance with claim 1, wherein said steps (4) through (6) 
are repeated one or more times until the "intervals smaller than said second predetermined 
intervals" in step (6) reach third predetermined intervals. 

4. (amended) A method in accordance with claim 1, wherein said predetermined interval 
is determined by using the following equation. 

Said first predetermined intervals are j^'^dv. Here, ; is a single-digit natural 
number, rfv is the increment in wave number, and kmax is the largest natural number 
satisfying the relationship jf»^*2pdv ^ Vmax, where Vmax represents the maximum 
calculation range from the center of the absorption line, and p is a natural number for 
controlling the calculation precision. 

5. (amended) A method in accordance with claim 3, wherein said predetermined interval 

is determined using the following equation. 

The third predetermined intervals are j'^'"dv. Here, / is a single-digit natural 
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number, dv is the increment of the wave number, and kmin is the maximum non-negative 
decimal number satisfying the relationship j'^'"^"pdv <a{a being approximately y/4) where y 
is an approximate value of the full-width at half-maximum of the absorption line, and p is a 
natural number for controlling the calculation precision. 

6. (amended) A method in accordance with claim 1, wherein said second predetermined 
interval is determined by using the following equation. 

For the sub-function with the (k - kmin + l)-th smallest width, the predetermined 
interval is j^dv. Here, / is a single-digit natural number, dv is the increase in wave number 
and k is such that kmin ^k< kmax, kmax is a natural number satisfying the relationship 
jkmax^ip^^ < Ymax, where Vmax represents the maximum calculation range from the center of 
the absorption line, and p is a natural number for controlling the calculation precision. 
kmin is the maximum non-negative decimal number satisfying the relationship j'^'ydv <a{a 
being approximately y/4), where y is an approximate value of the full-width at 
half-maximum of the absorption line, and p is a natural number for controlling the 
calculation precision. 

7. (amended) A method in accordance with any one of claims 1-3, wherein said 
interpolation is calculated by quartering the interpolation interval, with the function values 
represented by ya, yb, yc and the derivative values represented by ya, yh' , yc at the points 
where the interpolation interval {xo, xi) is divided, using the function values yo, yi and 
function derivative values yo', yi' at xo, xi, with the below-given Equation (1) as a function 
value interpolation equation. Equation (2) as a function derivative value interpolation 
equation and e as a non-negative decimal fraction. 
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8. (amended) A method in accordance with any one of claims 1-3, wherein said 
interpolation is calculated by dividing the interpolation interval into five parts, with the 
function values represented by ya, yb, yc, yd and the derivative values represented by ya, yb', 
ye, y/ at the points where the interpolation interval {xo, xi) is divided, using the function 
values yo, yi and function derivative values yo', y\ at xo, xi, with the below-given Equation (3) 
as a function value interpolation equation. Equation (4) as a function derivative value 
interpolation equation and £ as a non-negative decimal fraction. 
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9. (amended) A method in accordance with any one of claims 1-8, for increasing the 
speed of line-by-line calculations, when assuming the Voigt function to be K(x, y) and the 
shape of the absorption line displaced from the Voigt profile of the absorption line to be K(x, 
y) +/(^)/ by replacing 
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K{x,y) 

with: 

K{x,y)^AK{x,y) + Bf{x) 

and 

dK(x,y) 
dx 

with 

dK{x,y) ^^ dK{x,y) ^^dfjx) 
dx dx dx 

10. (amended) A method in accordance with any one of claims 1-8, for increasing the 
speed of line-by-line calculations, when assuming the Voigt function to be K{x, y) and the 
shape of the absorption line displaced from the Voigt profile to be K{x, y)f(x), by replacing 

K{x,y) 
with: 

K{x,y) = K{x,y)f{x) 

and 

dK{x,y) 
dx 

with 

^l^='-^f{,).K{x,y)^ 
dx dx ^ ^ dx 

11. (amended) A method in claim 10, for increasing the speed of line-by-line calculations, 
by using, in a sub-Lorentzian correction, 

Kix^y) = K{x,y)A exp(- B\x\) 
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and 



^ ^ ^1^)^ ^^^^^^^ J] 

12. (amended) A method in accordance with claim 9, for increasing the speed of 
line-by-line calculations, by replacing, in a line-mixing correction, 

with: 

K{x,y) = AK{x,y)+ BL{x,y) 

and 

dK{x,y) 
etc 



with 



dK{x,y) ^ 
dx 



{Ax + By)K{x,y)- {Ay - Bx)l{x, y)-^ 

Here, L(x, y) is the imaginary component of the function ti;(z) (the real part is ttie Voigt 
function), where complex number z = x + iy, defined by the following equation. 

w(z) = ^ £^^t£)j, = exp(-z')erfc{-iz) = K{x,y) + iL{x,y) 

(erfc{z) is a complex complementary error function). 

13. (amended) A program for increasing the speed of line-by-line calculations of multiple 

overlapping absorption lines, said program comprising: 

(1) a step of dividing a domain of a Voigt function representing the shape of an 
absorption line into a first range around the peak of the Voigt function and a skirt 
portion not contained in the first range, replacing the first range with a cubic function 
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whose function values and derivative values match with those of said Voigt function at 
the connection points with said Voigt function, and calculating the values and 
derivative values of said cubic function and the Voigt function in the skirt portion for 
each of first predetermined intervals; 

(2) a step of adding together the results of step (1) for a plurality of absorption lines; 

(3) a step of interpolating the results of step (2) to calculate the function values and 
derivative values in second predetermined intervals smaller than said first 
predetermined intervals; 

(4) a step of dividing said first range into a second range near the peak and a skirt 
portion not contained in the second range, replacing said second range of a "function 
representing the difference between the Voigt function and said cubic function" with a 
second cubic function whose function values and derivative values match with those 
of said "function representing the difference between the Voigt function and said cubic 
function" at the connection points with said "function representing the difference 
between the Voigt function and said cubic function", and calculating the values and 
derivative values of said second cubic function replacing said second range and the 
"function representing the difference between the Voigt function and said second 
cubic function" in the skirt portion for said second predetermined intervals; 

(5) a step of adding the results of step (4) to the results of step (3) for a plurality of 
absorption lines; 

(6) a step of interpolating the results of step (5) to calculate the function values and 
derivative values for intervals smaller than said second predetermined intervals; and 

(7) a step of adding the values of the "function representing the difference between 

the Voigt function and said second cubic function" to the results of step (6) for a 

plurality of absorption lines in said second range. 
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14. (amended) A program in accordance with claim 13, wherein the step of calculating 
function values and derivative values by interpolation is repeated while narrowing the 
interval until it becomes an interval of a minimum unit. 

15. (amended) A program in accordance with claim 13, wherein said steps (4) through 
(6) are repeated one or more times until the "intervals smaller than said second 
predetermined intervals" in step (6) reach third predetermined intervals. 
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